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Abstract 

The magnetic properties of a system of coexisting localized spins and conduction electrons are 
investigated within an extended version of the one dimensional Kondo lattice model in which 
effects stemming from the electron-lattice and on-site Coulomb interactions are explicitly included. 
After bosonizing the conduction electrons, is it observed that intrinsic inhomogeneities with the 
statistical scaling properties of a Griffiths phase appear, and determine the spin structure of the 
localized impurities. The appearance of the inhomogeneities is enhanced by appropriate phonons 
and acts destructively on the spin ordering. The inhomogeneities appear on well defined length 
scales and can be compared to the formation of intrinsic mesoscopic metastable patterns which are 
found in two-fluid phenomenologies. 

PACS numbers: PACS No. 71.27.+a, 71.28.+d, 75.20.Hr 
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I. INTRODUCTION 



The interplay of spin, charge and lattice degrees of freedom has been investigated inten- 
sively in many transition metal oxides and especially in perovskite manganites, which have 
recently attracted new interest due to the discovery of colossal magnetoresistance (CMR). 
The initial understanding of the properties of manganites was based on the double-exchange 
mechanism within a Kondo lattice.^ However, neutron scattering^ and electron diffraction^ 
experiments have revealed the simultaneous presence of charge and spin superstructures 
which makes these early theoretical approaches incomplete. 

The existence of charge, lattice and spin modulations with a doping-dependent wave vec- 
tor in La2_xSrj.Ni04,i^ or the presence of charge ordering at half filling in Ndo.sSro.sMnOg^ 
and similar compounds, such as Pro.sCao.sMnOs,"^'^'^ all suggest that the theoretical un- 
derstanding has to be extended to account for effects stemming from the lattice in order 
to understand the doping dependent phase diagram and the richness of phases that are ob- 
tained. More recently, in electron doped charge ordered manganites Lai-i^CaajMnOs charge, 
orbital, and magnetic ordering has been observed- for the first time. Several experiments 
also confirmed this behavior in other compounds . Charge ordering has also been found 
for other doping systems as in Bii_a.Caa;MnOar^ and in Lai_i.Ca2^Mn03 (doped with P^)^ 
for X > 1/2. 

Because of these complex behaviors, manganese oxides have been intensively studied 
during the last few years. Since realistic models of these perovskites is impossible to solve 
completely, different approximation schemes have been introduced to account for the indi- 
vidual properties of their rich phase diagram. Typically, the electronic degrees of freedom 
are described by a Kondo lattice model (KLM) which in the strong Hund coupling limit 
reduces to the Zener— double-exchange hamiltonian. Early theories- proposed to explain 
the physics of CMR materials have focused primarily on this model. As the previously noted 
experimental findings show, all the approaches based entirely on the phenomenon of double- 
exchange are incompleted^. Understanding the complex phase diagram of the perovskites 
can be resolved by including additional physics in the Kondo lattice, in the form of electron- 
phonon coupling^- originating from a Jahn- Teller splitting of the Mn^+ ion. Consequently, 
our first goal here is to understand the lattice effects in the KLM by including explicitly the 
interaction with the lattice degrees of freedom. 
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A further consequence of the strong Hund couphng is that three locahzed orbitals of 
the Mn^+ ions will be aligned, giving rise to a 5* = 3/2 localized spin. As most of the recent 
approaches to CMR are based on Monte Carlo simulations of one dimensional models,— 
it has been argued that for finite temperatures a classical spin represents a reasonable ap- 
proximation. However, this may not be the best approach, as the quantum fiuctuations are 
the strongest in one dimension. Hence, our second goal here is to develop a bosonization 
approach which solves the full quantum spin KLM in one dimension. For completeness we 
will allow for both ferro- and antiferromagnetic coupling in the KLM. In CMR materials 
where the coupling is ferromagnetic, it is obvious that double-exchange dominates. We will 
show in Section II that this is also the case for antiferromagnetic coupling, as in, e.g., heavy 
fermion compounds^^. 

The paper is organised as follows, in Section II we investigate the presence of double- 
exchange ferromagnetism in both J < and J > cases. Section HI contains a description of 
our starting hamiltonian. Section IV contains the details of the bosonized solution. Section 
V is devoted to a comprehensive description of the localized spin ordering. In Section VI 
the obtained phase diagram is analysed in detail. Section VII contains the conclusions. 

II. DOUBLE-EXCHANGE 

In order to gain a more transparent understanding of the double-exchange interaction we 
investigate first the case of two sites and one conduction electron, as was done originally by 
Anderson and Hasegawa, and de Gennes.— In the case of ferromagnetic coupling (J < 0) 
the ground state energy is -Eq, j<o = —\J\/^ " ^ with wave function |'?/'o)j<o = |'^de)j<o, ^ = 
I fTzTz! fTz 0) + I 0, fTztz)) where fl^^ and ]z refers to the z component of the impurity 
and conduction electron spins, respectively. As can be seen, ferromagnetism arises here via 
an Ising type coupling, which allows for description of the ground state within a simple 
semiclassical approximation.i^ 

For J > the situation is completely changed due to the singlet formation of localized 
and conduction electron spins. Due to this Kondo singlet formation, which is absent for the 
J < case, double-exchange is often ignored in the discussions of the J > KLM. This case 
has usually been discussed in terms of the competition between Kondo singlet formation and 
the RKKY interaction. For a half filled band, this point of view appears to be sufficient. 
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However for a partially filled conduction band an overwhelming amount of numerical data^i 
proves that ferromagnetism appears for stronger coupling. 

This cannot be explained in terms of RKKY, which operates at weak-coupling, nor in 
terms of Kondo singlets, since they are non-magnetic. The missing element is double- 
exchange orderingiSi due to an excess of localized spins over conduction electrons. Double- 
exchange requires only that the number of conduction electrons is less than the number of 
localized spins. It operates in any dimension, for any sign of the coupling J, and for any 
magnitude Sj of the localized spins.—. The double-exchange interaction is specific to the 
Kondo lattice, and is absent in single- or dilute-impurity systems in which the situation is 
reversed, and the electrons greatly outnumber the localized spins. 

Double-exchange is conceptually a very simple interaction: each electron has on average 
more than one localized spin to screen, and consequently hops between several adjacent 
spins gaining screening energy at each site, together with a gain in kinetic energy. Since 
hopping is energetically most favourable for electrons which preserve their spin as they hop 
(called coherent hopping), this tends to align the underlying localized spins.-^. 

For two sites, this causes a mixing of the total spin and an enhancement of the Hilbert 
space, where now 16 elements have to be considered. The ground state energy is given 
by -Eo, j>o = - J/4: - + 2Jt + 4t^/2 with wave functions \ipo)j>o oc \ipKs)z + - 
^o,j>o)]{| ^4.,^. 0) + | 0, itdz)-\i[zhAz 0)-| ^,0,^4,)}, where the Kondo singlet 
iV'Ks)^ states are | il^l^, 0) - | Jj-J^, 0) + | 0, il^lz) - I 0, J|4^). |^/'o)j>o involves 
six basis elements (the degeneracy is partially lifted by conduction electron hopping) and 
hence falls outside the four dimensional space needed to establish double-exchange for J < 0. 

In order to invoke double-exchange as well, all three spin directions, x, y, and z, have to 
be considered: |V'o)j>o oc [1 - l/(J/4 - ^o, j>o)] \^ks)z + [l/(<//4 - ^o, j>o)] U^be) .j>o, x + 
|V'de)j>o,s/ + |^de)j>o,^}, where IV'de) J>0, a=x or y 

, 4a 0) + I ^„ 0, \^a]a)}/V2 and |^DE)j>o,a=. = | tziz, tz 0) + \ tz 0, ^4.)- In spite of 
this extra complication, it is apparent from the above that also for J > 0, ferromagnetism 
appears. 

We show in Section V that coherent conduction electron hopping over a characteristic 
length A may be incorporated into a bosonization description which keeps the electrons 
finitely delocalized. At lengths beyond A, the electrons are described by collective density 
fluctuations, common to one dimensional Fermi systems. The electrons remain finitely 
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delocalized over shorter lengths, and describe coherent hopping over several adjacent sites. 
This tends to align the underlying localized spins at stronger coupling. A measures the 
effective range of the double-exchange interaction, and hence is proportional to the width 
of the magnetic polarons. 

III. THE MODEL 

The KLM considers the coupling between a half-filled narrow band (localized d or /) 
and conduction electrons. Even though studied intensively for the last two decades, the 
understanding of the KLM remains incomplete. Only in one dimension have numerical 
simulationsii and bosonization technique a^^i^^ been carried through to admit predictions 
about the phase diagram of the KLM. No investigations exist for the case where the KLM 
is extended to account for contributions stemming from the phonons, which is of special 
relevance to OMR materials. In particular, the small doping regime of these systems, which 
are ferromagnetic at low temperatures, seems to be appropriate for modelling within the 
KLM argumented by interactions with the lattice. In the following we present bosonized 
solutions of the KLM where on-site Coulomb and specific phonon contributions are explicitly 
included. This "extended" KLM model allows spin-and magnetoelastic-polaron formation, 
which we believe are of major importance in understanding these complex materials. 

The Hamiltonian of the KLM in the presence of on-site Coulomb interaction reads: 

^^KLM+U = -iim>S+l,<x + h-C.) 

+ JY. Sd,i-Scj + UY^rij^^rii^^ , (1) 

3 3 

where t > is the conduction electron hopping integral, Sdj = \Y.cr,cr' '^d,j,a^^,(T'Cd.j,a'j 
Scj = \J2a,a' c]^cr^<Ty(^j,a' ^ Pauh splu matrlccs. Fermi operators Cj^^,Cj^^ 

with subscript d refer to localized d-spins, while those not indexed refer to the conduction 
electrons. The on-site Coulomb repulsion is given by the Hubbard term proportional to U. 
In the CMR materials the localized states are represented by the threefold degenerate Mn 
(i-electrons with total spin 3/2. However, for reasons of transparency, the localized spin 
is approximated here by spin 1/2. The properties of the model are qualitatively independent 
of the magnitude of the localized spins, i.e., the basic features of the phase diagram for the 
three dimensional case appear in one dimensions as wel&^. 
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In the following the Kondo coupling J is measured in units of the hopping t and both 
cases, antiferromagnetic (J > 0) and ferromagnetic (J < 0) couplings, will be considered. 
The conduction band filling is given by n = N^/N < 1, where is the number of lattice sites 
and Nc is the number of conduction electrons. In order to understand the broad properties 
of CMR materials, we also allow for the number of impurity spins, A^^^, to vary, in such 
a way that Nd/N < 1. In this way, doped or dilute Kondo lattice systems can also be 
studied. It will be shown that a dilute Kondo lattice system is dominated by short range 
antiferromagnetic correlations. 

In one dimension, the electron-phonon coupling can either be of inter-site (Su-Schrieffer- 
Heeger (SSH)-^) or on-site (Holstein^^) character. The models we study assumes a disper- 
sionless phonon mode with frequency uj. The neglect of the dispersion of bare phonons is 
not essential since it is absent in the Holstein model and the acoustic phonons are decoupled 
from the low energy electronic spectrum in the continuum limit of the SSH model.— In this 
approximation the Holstein coupling to dispersionless phonons has the following form: 

i^Holstein = {Pli^ci + + ' (2) 

i 

with the conduction electron density n^^ at site i, the lattice displacement q^, its conjugate 
momentum p^, the electron-phonon coupling strength /?, the spring constant K and the ionic 
mass M. Within the SSH model the coupling to phonons is modified to 

where denotes the electron-phonon coupling strength. Thus, the starting hamiltonian is 

H = HxLM+V + -f^Holstein + -f^SSH • (4) 

These phononic contributions may not describe the full complexity of the phononic cou- 
plings observed in real materials because of the phase space constraint of any one dimensional 
calculation. However, the results capture the essence of the Kondo lattice coexisting with 
phonons, and being an exact solution, it represents a vital source of information due to the 
lack of comparable solutions applicable to colossal magnetoresistance materials. 
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IV. THE EFFECTIVE HAMILTONIAN 



A large class of one dimensional many-electron systems may be described using bosoniza- 
tion techniques'^: The electron fields may be represented in terms of collective density 
operators which satisfy bosonic commutation relations. Bosonic representations provide a 
non-perturbative description which, in general, by far easier to evaluate than a formulation 
in terms of fermionic operators. 

The underlying bosonization scheme follows standard procedures'^ by first decomposing 
the on-site operators into Dirac fields, c^. ^ = J2t e^'^^^'^r,a{x), where kp = 7m/2, with spinor 
components r = ± (+/- being the right/left movers) and kp = 7m/2. Next we bosonize the 
Dirac fields with '^r,a = exp(i$^^o-) / VSvrA, where 1/A is the ultraviolet cutoff. For the scalar 
Bose fields, ^r,a{x), and their conjugate momenta, Ilr,a{x), ^r,a{x) = J!!^^dx'Ilr,a{x'), are 
used in standard Mandelstam representation by means of which a momentum cutoff via the 
Fourier transform is introduced A(fc) = exp(— A|fc|/2). If the distance between the impurity 
spins is larger than A, the electrons will behave as collective density fluctuations.— Thus, the 
Fermi fields can be represented in terms of density operators which satisfy Bose commutation 
relations: 

c^,^^^ = expliTkpx) expi{6p{x) + T(j)p{x) + (r[6^{x) + T(f)^{x)]}/2 , (5) 
where the Bose fields for u = p,a are defined by 

Mx) = i{ii/N) J2 e*'"K(A:) ± iy.{k)]A{k)/k , (6) 

with + corresponding to the number fields ipu = <Pu and — to the current fields ipi, = 9^. 
The charge (holon) and spin (spinon) number fiuctuations are defined as Pr{k) = J2a PT,a{.k), 
and (Jr{k) = Z^o- ^Pt,(t(^)- All rapidly oscillating terms originating from e.g. backscattering 
and umklapp processes are neglected, since they contribute only at exactly half filling. The 
localized d electrons can neither be bosonized nor Jordan- Wigner transformed since no direct 
interaction exists between them. 

Substituting these representations into Eq. gives the bosonized KLM hamiltonian 

+ A E {^os[Mj)] + cos[2kpj + 0,(j)]} (e-^'-(^')5+ + h.c.) 
j 
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^ Y: sin[0.(j)] sin[2A;^j + 0p(j)]^d„- • (7) 



The charge and spin velocities are given below. 

Considering the phononic contributions, in standard bosonization language, Eq. El sim- 
plifies to + iff -P'' + Hf-^^, where 



p 

+ ydx[ulix)+uXix)] , (8) 

and uq = = ^K/M are their respective phonon frequencies. The electron-phonon 
forward scattering term is simply 

^2^-"" = 72^ E[P+(-p) + P-{-p)]Mp) ■ (9) 

iV p 

On the other hand, the rapidly oscillating phonon-assisted backward scattering term will 
acquire an extra factor exp[±i7TSnx] = exp[±i{2kF — 7r)x], in the form: 

Hf'' = 71 E / dx[^l^^_^^^e^-'^'n^^{x) , (10) 

with 7i = 72 = /5/a/M, where we used the same subscripts for backward and forward 
scattering as in g-ology.-^ 

In the continuum limit, the SSH term in contrast to the Holstein coupling, gives only 
two terms H^^ + _ the standard phononic component and a rapidly oscillating back 

scattering term. H^^ is given in Eq. |H1 while the back scattering term differs from 

Eq. ^Jonly in a form factor: 

i/l'i-P*^ = 7-1 E / rfx[^z/vl/t^^vl;__^e-^"1$.(x) , (11) 

with 7_i = Aacr/ v^M. However, the fact that the forward scattering term is missing means 
that the SSH coupling will not give any contribution to the effective hamiltonian away from 
half filling. 

Thus, the transformed hamiltonian of Eq. (0} is: 

^ = 2^E {^lU) + [d.Mj)?} + ^E[n^(p) + ^o*o(p) 
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+ ^ E {cos[0.(j)] + cos[2A;^j + (e-^'^^^^^J^- + h.c.) 

- 4^ E MM)] sin[2fcFj + <l^mSh + ^ E[^x0.(j)]'^d,, • (12) 

J 3 

If holes are present in the array of d-spins, all terms proportional to S are zero. The charge 
and spin velocities are 



Vp 



= vp^Jl + U/t:vf - P^/ttKvf , = vp^l - U/nvp + P^/nKvp , (13) 



where the Fermi velocity is vp — 2sin(7rn/2) in units of t. 

It is important to note that a renormalization of the spinon-holon velocities appears 
here due to the Hubbard and phonon terms which act oppositely on the corresponding 
velocities. While the Hubbard term leads to a localization of the spinous and an increased 
hopping of the holons, thus supporting a magnetic ground state, the phonons delocalize the 
spins, but localize the charges and act destructively on the magnetic properties. It is worth 
mentioning that the Hubbard term alone already suffices to establish two time scales for the 
holon-spinon dynamics. But an important renormalization of the critical properties of the 
system is achieved through the variable phonon coupling, which, as will be shown below, 
estabhshes the existence of a Griffiths phase. The competition between the Hubbard and 
the phonon term obviously vanishes for U = j3'^ / K. 

In the following, effects arising from the localized spin d impurities, double-exchange, 
the phonons and Hubbard interactions will be discussed in more detail. The localized spin 
d impurities act via double-exchange on the hopping electrons so as to preserve their spin 
when moving through the lattice in order to screen the localized spins which are in excess of 
the conduction electrons, i.e. N < > N^. This, in turn, leads to a tendency to ahgn the 
localized spins and results in an additional screening energy for the conduction electrons. 



V. LOCALIZED SPIN ORDERING 



In order to determine rigorously the phase diagram and investigate the ordering of the 
local spins due to the formation of polarons, we firstly apply a unitary transformation. This 
is the simplest method for determining the ordering of the localized spins induced by the 
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conduction electrons and it has been used successfully elsewlierei^>i^. This is achieved by 
choosing a basis of states for the unitary transformation in which competing effects become 
more transparent, i.e. a transformation which changes to a basis of states in which the 
conduction electron spin degrees of freedom are directly coupled to the localized spins. 



to Eq. ()12j) up to infinite order, thus avoiding truncation errors. 

Secondly, we explicitly take into account the Luttinger hquid character of the Bose fields, 
i.e., use their non-interacting expectation values such that the effective Hamiltonian for the 
local spins is derived as: 



the site j and subtracts all those to the left of j: K{j) = {J/2vf) Z]«Si('S'dj+/ — S^j_i). This 
term gives the crucial difference between the Kondo lattice and dilute Kondo lattice, as will 
be explained in the following. 

The most important term in Eq. ()14|) is the first one. 



This term shows that a double-exchange ferromagnetic interaction appears for both J > 
and J < coupling and even in the dilute Kondo lattice model. This coupling is non- 
negligible for Nd > Nc and j — j' < A and its strength decreases with increasing distance 
between impurity spins. The value of JgfT is plotted for different values of U and /5 in Fig. [H 

Fig. ^ shows clearly that the Hubbard interaction enhances double-exchange and conse- 
quently Jeflf. As noted earlier, see the discussion following Eq. (fTSj) . this is due to the fact 
that the Hubbard term leads to a localization of spinous and an increase in the hopping 
of holons, hence increasing the double-exchange. On the other hand, the electron-phonon 
interaction counteracts this effect by localizing the holons and thus decreasing Jes- 

Fig. ^ also shows that the length A characterizes the effective range of the double-exchange 
interaction. This interaction originates from the bosonization of the conduction electron 



Correspondingly, we choose the operator S = i{J/2'K)^JvF/v^ J2j ^cr(j) S^j, which is applied 
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FIG. 1: The range of the ferromagnetic interaction Eq. (|15|) for different values of U and (3. JeS 
is the interaction strength in units of j'^v'^/Air'^vp 



band in the following way: at wavelengths larger than A, the electrons combine to form 
collective density fluctuations, which involve large numbers of electrons satisfying bosonic 
commutation relations. This is the standard behavior of one dimensional many-electron 
systems for weak interactions^^. At wavelengths smaller than A, the density fluctuations are 
not collective, and loose the bosonic character. Since bosonization only applies to fluctua- 
tions beyond A, the bosonization description is equivalent to keeping the electrons finitely 
delocalized within the range of A, with the electrons preserving their spin over this range. 
Eq. (fT3j) describe the ordering induced on the localized spins by the finitely delocalized elec- 
trons. Thus, A corresponds to the effective delocalization length related to the spatial extent 
of the polarons (polaron radius), i.e., the effective range of double-exchange, as shown in 
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FIG. 2: A snapshot of the magnetic polarons. (4) and | (|) refers to the z component of 
the impurity and conduction electrons spins. The dashed curve represents the spin domain walls 
(kink-antikink pairs). The electrons will oscillate between these domains walls (dotted circle) and 
hence couple the impurity spins through double-exchange (for details see text). A is the average 
length of the polarons, i.e., the effective range of the double-exchange. 



Fig. El 

Double-exchange becomes inefficient if the distance between the impurity spins is larger 
than A. In general A has a complex behavior, depending on J, [/, /? and N^- For large 
distances and low density, Eq. (jl5p gives an effective range for the double-exchange which 
decays exponentially with a characteristics scale y^2/J. Hence, in the following we will use 
for A its low density (e.g., in the example depicted in Fig. 12)) value: A ~ y^2/J. And since the 
interaction Eq. ()15|) is short range for all finite A we approximate it by its nearest-neighbor 
form: 

j^jj = —J^ / dkcoskA'^(k) . (16) 
Jo 

Here we emphasize again that Eq. (jlfij) is valid for both J > and J < couplings. 
Since the J > case is usually not considered in discussions of the KLM, we emphasize the 
following points: (1) the term originates, via bosonization and then the unitary transforma- 
tion, from the kinetic energy term —t I]j>(Cj,o-Cj+i + h.c.) and the forward scattering part 
of ( J/2) X^jl'^jT ~ (%o- ~ ^]crCja) ill the KLM hamiltonian Eq. ((T)). (Note that the 

Bose representations for the electrons in these terms are exact.) (2) Eq. is independent 
of the sign of J, and takes the same form for any magnitude Sj of the localized spins. (3) 
Since Eq. p6|) is of order J^, whereas the remaining terms in the transformed hamiltonian 
Eq. p4|) are of order J, the interaction Eq. (fT^ dominates the ordering of the localized spins 
as J increases. All these properties are identical to those of a double-exchange interaction. 
This leads us to identify Eq. (|T6|l as the double-exchange interaction in the KLM. We also 
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emphasize that Eq. p6p is ferromagnetic for all (possible) choices of the cut-off function 
A{k). 

Based on property (1) above, a simple characterization can be given for the double- 
exchange valid at low conduction band filling. The simple hamiltonian term (1) from above 
satisfy a standard nonlinear Schrodinger equation: dl^pai^) + {Jmci/2)\ipf^{x)\'^4'a-{x) = 
2meiEipcr{x) {rriei. being the bare electron mass, and ilJrrix) the electronic wave function) 
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with soliton solutions ^jjaix) oc e*^sech(a;y Jmei./4); see Ref. 

These soliton solutions correspond to spin domain walls of finite size (kink-antikink pairs) 
and lead to a gain in electronic energy of —a for antiferromagnetic coupling, and of +cr for 
the ferromagnetic shown in Fig. |21 Physically the solutions resemble the dressing of 

the electron by a finite range of parallel (antiparallel) local spins and consequently can be 
identified with polaronic type objects. From the previous considerations it can also be con- 
cluded that, when including the interactions with the phonons, the tendency towards charge 
localization is enhanced and increases this polaronic effect. Since the lattice also experi- 
ences a renormalization due to the coupling to the electronic degrees of freedom, substantial 
ionic displacement patterns will develop and the formation of magnetoelastic polarons takes 
place. Similar results are obtained by decoupling electronic and phononic degrees of free- 
dom through a homogeneous Lang-Firsov transformation, where the localization stems from 
band narrowing. In accordance with previous results, the polaron radius is characterized by 



a length scale proportional to y 2/ J . This new length scale differs from the free conduction 
electrons mean free path and gives rise to competing time scales: slow motion of the pola- 
ronic carriers and fast motion of the free electrons, thus providing dynamics of two types 
of particles and a close analogy to a two fluid scenario.^ Since the polarons are in general 
randomly distributed within the local spin array, these states can be viewed as intrinsic 
inhomogeneities involving spin fluctuations and short-range spin correlations. In addition, 
these new slow dynamics will exhibit a peak in the spin structure factor at 2kF — vr instead 
of the simple 2kp RKKY signal. A similar observation has also been made^ using numerical 
approaches. 
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VI. THE PHASE DIAGRAM 



For the Kondo lattice model, K{j) as the number of localized spins to the left and 
the right of a given site j are the same. If, however, we have a small concentration of holes 
in the array of localized spins, then - opposite to the previous case - K{j) is non-vanishing 
since the hole spins are no longer necessarily equally distributed to the left and the right 
of a given site. This yields K{j) ^ {—iy{J/2vp), which gives rise to a staggered field and 
antiferromagnetic ordering.— 

Since our main interest here is to explore the occurrence of ferromagnetism in the presence 
of the Hubbard and phonon terms, we take K{j) vanishingly small and focus on the transition 
between the paramagnetic and the ferromagnetic state controlled by the first two terms of 
Eq. (fT^ . Hence Hcs reduces to a quantum transverse- field Ising chain. This model without 
backscattering is known^° to undergo a quantum phase transition from a ferromagnetic to a 
paramagnetic phase. In the case of Eq. (fT^ . where a backscattering is also present, a similar 
phase transition occurs as the coupling J is decreased: the conduction electrons become less 
strongly bound to the localized spins, and tend to extend over spatial ranges beyond the 
effective range A of double-exchange ordering. Double-exchange becomes less effective, i.e., 
the magnetic polarons are loosely bound and regions of ordered localized spins begin to 
interfere as the conduction electrons become extended. The interference leads to spin-flip 
processes given by the transverse field in H^g. 

The transverse-field, h{j) = (J/27rA)[l + cos(2/i;i;'j)], in Hes, includes two low-energy spin- 
flip processes by means of which the conduction electrons disorder the localized spins. One 
spin-flip process is backscattering, and is accompanied by a momentum transfer of 2kp from 
the conduction electrons to the localized spins. Since the chain of localized spins will tend 
to order so as to reflect this transfer, the transverse-field corresponding to backscattering 
spin-flips is sinusoidally modulated by 2k p. The other low-energy spin-flip process in ifgff 
is forward scattering. This involves zero momentum transfer to the localized spins, and the 
corresponding transverse- field is a constant (i.e. has modulation zero). 

However, away from half filling the transverse field h{j) will have an incommensurate 
modulation 2kp with respect to the underlying lattice of localized spins. Hence, the conduc- 
tion band is unable to either totally order or totally disorder the lattice as the ferromagnetic 
to paramagnetic transition occurs. There remain dilute regions of double-exchange ordered 
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localized spins or magnetic polarons in the paramagnetic phase as only a quasi-commensurate 
fraction of the conduction electrons become weakly-bound, and become free to scatter along 
the chain. The remaining ordered regions are sufficiently dilute to prevent long range corre- 
lations, but their existence dominates the low-energy properties of the localized spins near 
the transition. 

These considerations motivated us to treat the transverse field as a random variable. The 
effective Hamiltonian can thus be replaced by a random transverse field Ising model:— 

j j 

where Jefr is given in Eq. ()16|) and the ferromagnetic coupling strictly vanishes if i > A. 
The random fields, hj, are generated by (1 + cos[2k pj]) at large distances, where cos[2k pj] 
oscillates unsystematically with respect to the lattice. The large values cos[2kpj] ^ 1 which 
are responsible for spin fiips, are then well separated and are driven by a cosine distribution 
analogous to spin-glasses.— 

Using the extensive real space renormalization group results of this model by Fisher, Ref. 
(311). we determine the location of the quantum critical line describing the paramagnetic to 
ferromagnetic transition at 

^crit = (7rV4)sin(7™/2) {1 - U/[2nsin{7in/2)] + (3y[2nK sm{7Tn/2)]y/^ . (18) 

For values J < Jcrit a paramagnetic state exists which is dominated by polaronic fiuctuations. 
For J > Jcrit ferromagnetism appears. The transition between these phases is of order- 
disorder type with variable critical exponent 6 = ln( Jcrit/ J). 

The behavior described by JTcrit is simply understood in terms of magnetic polarons 
introduced previously, see also Fig.|21 Considering for the moment the standard KLM model, 
a summary of the known results^ is presented in Fig. El Reducing J from intermediate 
values in the ferromagnetic phase, the infinite cluster characterizing strong ferromagnetism 
is broken up into several large clusters as the quantum fiuctuations hj, controlled by the spin- 
fiip interactions, become stronger. The individual clusters are the spin polarons depicted 
in Fig. 121 These magnetic polarons are weakly ordered in this phase, i.e., form a polaronic 
liquid which exists for —0.7 < 6 < hence with the upper boundary determined by J ~ 
2Jcrit, see Fig. |31 This is not a true transition line, but rather marks the crossover to a 
Griffiths phase^^ characterized by singularities in the free energy over the whole range of 
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FIG. 3: The ground-state phase diagram of the standard KLM. The sohd (critical) line is from 
Eq. H18|) . The dashed lines separate strongly ordered (conventional ferromagnet) / disordered 
(conventional paramagnet) phases from their weak Griffiths phase (polaronic liquid) counterparts. 

8. For small 8 the correlation length is ^ ~ beyond which the system is ordered. 
The spontaneous magnetization Mq oc \8\^ ■, with 7 = (3 — -\/5)/2 ^ 0.38,-^ while for small 
applied fields H the magnetization M(H) oc Mq[1 + 0{H'^^^^S In H)]; the susceptibility is 
infinite with a continuously variable exponent. The mean correlation function is given by 
(^/x)^/6e-^/«exp[-3(7rx/0^/^] for x > C-— 

Further lowering J, we reach the true phase transition Eq. p8|) . The correlation length is 
infinite, the magnetization M{H) oc | lnH\^"' for small H, and the mean correlation function 
is critical x"'''. The transition curve for different values of U and f3 is shown in Fig. |3 As 
argued previously, the Hubbard interaction makes the ferromagnetic phase more robust, as 
it increases the strength of the double-exchange and as such the length of the magnetic 
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FIG. 4: The quantum critical line Eq. describing the ferromagnetic transition for different 
values of U and /?. For each of these examples the ferromagnetic phase is above the phase transition 
curve. 

polarons. 

Calculating the effective length of polarons, A on the transition line, see Fig. supports 
this finding. The length of the polarons increases strongly with ?7, hence making the ferro- 
magnetic phase more difficult to break up. However, the consequence of this effect is that 
the Griffiths phase below the transition will diminish strongly or vanish in most cases. 

The effect of phonons is opposite to this. The phonons act destructively on the magnetic 
properties. The ferromagnetic phase becomes smaller, see Fig. |3] and the length of the 
polarons are slightly decreased in the presence of the phonons as shown in Fig.El But, as will 
be explained later, the boundaries of the lower Griffiths phase will remain the same. Thus, 
contrary the the effect of t/, where the Griffiths phase vanishes, for strong phononic couplings 
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FIG. 5: The effective range A of the double-exchange interaction in units of the lattice spacing 
against filling n on the critical line. 

the magnetoelastic polarons extend over a much larger phase space and will dominate the 
phase diagram of the KLM. 

Immediately below the critical line {6 > 0), the system represents a weakly disordered 
Griffiths phase. The remaining polarons occupy a small fraction of the system length but 
behave as if they are still in the ordered phase; their magnetization 6^ per unit length is 
identical to Mq of the weakly ordered phase.— These remaining rare polarons dominate 
the low-energy physics. Hence, crossing the phase transition line disorders the polarons. 
The transition is indeed an order-disorder transition of the magnetic polarons. This regime 
can be viewed as a paramagnet with locally ordered ferromagnetic regions. This scenario 
resembles a two-fluid picture, i.e., a polaronic liquid, with intrinsic inhomogeneities which 
involves spin fluctuations and short-range spin correlations. 

In this weakly disordered phase, the magnetization M(H) oc S^{H^^[5\n(l/H)+const.] + 
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FIG. 6: The ground-state phase diagram for U = and /3 = 5. The sohd line is true phase 
transition line from Eq. H18|) . The dashed lines separate the conventional ferromagnet / para- 
magnet phases from their weakly disordered (polaronic liquid) counterparts. Note the extent of 
the polaronic liquid phases compared to the U = and /? = standard KLM (shaded region) 
corresponding to Fig. Q. 

0[H*^6\n{l/ H)]}; thus M{H) has a power law singularity with a continuously variable 
exponent 27; as in the weakly ordered phase the susceptibility is infinite. 

The mean correlation function decays less rapidly than in the weakly ordered phase, but 
takes the same form (^/x)^/6e-^/«exp[-3/2(7rx/0^/^] for a; > ^ = 1/6^. According to ifcrit, 
the weakly disordered Griffiths phase extends down to J = 0. However, as the disorder 
increases, the third term in Hes is no longer negligible. At very low J, the last two terms 
in will dominate; this corresponds to free spins in a field with dominant correlations at 
2/cj? of the conduction band, and is responsible for the observed peak in the localized spin 
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structure factor.^- This is the strongly disordered conventional paramagnetic phase, see Fig. 

m 

As the boundary of this Griffiths phase is governed by the full i^efr from Eq. (jl4p rather 
than Jcrit only, it will not be effected by U and (3. This explains why the extent of this 
Griffiths phase increases dramatically in the presence of phonons and vanishes for large U. 
Thus the phonons indeed will enhance charge localization and increase the polaronic effect 
and act destructively on the spin ordering. In the presence of phonons, these magnetoelastic 
polarons will dominate a larger phase space below and above the ferromagnetic transition 
compared to the standard KLM. This is shown explicitly in Fig. where the (3 = 5 case is 
compared to the (3 = 0, and [/ = in both cases. It can be seen that the polaronic liquid 
occupies a much larger phase space in the ffist case compared to the standard KLM (shaded 
region). This proves that, at least in the one dimensional KLM the local inhomogeneities 
are enhanced by phonons, creating magnetoelastic polarons similar to intrinsic mesoscopic 
patterns from two-fluid model phenomenologies. 

It is interest ing to mention that similar behavior is expected in higher dimensions also. As 
shown in Ref. (1341 ). the randomness generated through local inhomogeneities is the driving 
force, rather then dimensionality, in any real space renormalization approach. Thus most 
of the properties presented previously will survive in higher dimensions. This have been 
confirmed also by extensive numerical calculations.— 

An important prediction of our theory is the finite temperature behavior of the polaronic 
liquid phase below the transition hne, the so-called weakly disordered paramagnetic phase. 
This is the region of the phase diagram which is enhanced by phonons and as such dominated 



by magnetoelastic polarons, as explained previously. Here, following closely Ref. ()3lD. we can 
determine the finite temperature susceptibility. At small temperatures the weakly coupled 
polarons behave as free spins with momenta 6'^~^ \n{T*/T) and the susceptibility will become 
Xpoiaron ~ 5^''' (T/T*)^*^"^ lu^ (T/T*) . T* IS a temperature scale at which the renormalization 
flow is stopped when the temperature fluctuation scale equals the quantum fluctuation 

The susceptibility diverges at T in the Griffiths phase as expectecP^. However, if 
6 increases, e.g., 5^1/2 the system moves away from criticality, the susceptibility at zero 
temperature will be finite and will no longer be dominated by the large rare polarons but 
rather by the more typical smaller ones, becoming a strongly disordered phase. On the 
temperature scale, this corresponds to increasing the temperature. Thus, at higher temper- 
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atures the susceptibility will become Xpoiamn ~ {T/T* - 1) ln(T/T*) = (1 -T/T*) \n{T* /T). 
Interestingly, this is identical to the phenomenological form of the susceptibility that has 
been associated with heavy quasiparticles in all heavy fermion compounds.-- For these com- 
pounds T* is the temperature at which the heavy quasiparticle liquid start to emerge from 
a standard Kondo lattice behavior. Note that in our approach of deriving Xpoiaron "we moved 
away from the critical point, in agreement with recent calculations^^ that the heavy fermion 
systems are in the vicinity off, but not on a Griffiths singularity. 

VII. CONCLUSIONS 

In summary, we have derived an effective Hamiltonian from a one-dimensional Kondo 
lattice model extended to include effects stemming from coupling to the lattice and in the 
presence of an onsite Hubbard term, which accounts for the conduction electron Coulomb 
repulsion. The results are: i) A ferromagnetic phase appears at intermediate |J| due to 
forward scattering by delocalized conduction electrons, ii) Ferromagnetism is favoured by the 
Hubbard term, while it is suppressed by the electron-phonon coupling, iii) The paramagnetic 
phase is characterized by the coexistence of polaronic regimes with intrinsic ferromagnetic 
order and ordinary conduction electrons, iv) In the paramagnetic phase, two time scales 
compete with each other - reminiscent of a two-fluid model - and the variability of the critical 
exponents suggests the existence of a Griffiths phase. 

The results are related to the small-doping regime of CMR materials which are ferromag- 
nets at low temperatures, since here the coupling to the phonons has been shown to dominate 
the paramagnetic-ferromagnetic phase transition. Note that moving away from this criti- 
cality into the weakly disordered polaronic liquid phase the results, in particular the finite 
temperature susceptibility, show behavior characteristic to all heavy fermion compounds.— 

In regard to the CMR materials, it is interesting to point out the discrepancy between 
infinite dimensional calculations and the present one dimensional results. Many approxi- 
mate calculations to model CMR— have been made in dynamical mean-field theory, which 
is an infinite dimensional approximation and therefore incapable of capturing spatial in- 
homogeneities. In the present work we have approached the CMR materials via a one 
dimensional approximation, but with techniques able to describe fiuctuations of short-range 
order. Our results show that strong intrinsic spatial inhomogeneities of Griffiths type dom- 
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inate the behaviour of the Kondo lattice. Consequently the inhomogeneities exhibit clear 
statistical scaling properties as a function of the proximity to a quantum (order-disorder) 
critical point. The phonons enhance the inhomogeneities, which, in a good approximation 
behave as a supercritical (metastable) phase of a two fluid model. 

Even though various bosonization schemes have been used for the one-dimensional 
KLM ,^^i^^ non of the previous approaches took into account phonons. The inclusion of 
phonon degrees of freedom has been shown here to be relevant in creating local magnetic 
inhomogeneities. It is important to emphasize that the properties of the system are con- 
trolled by intrinsic inhomogeneities. This means that, in a renormalization group approach, 
the dimensionality should not matter .i^^^ Thus, a similar behaviour is expected in realistic 
two- and three-dimensions, which clearly merit further detailed study. 
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